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1*.  SPONSOPINO  MILITARY  ACTIVITY 

Office  of  Nava!  Research 

Department  of  the  Navy 

Arlington,  Virginia  22217 

An  analysis  based  on  the  lifting  surface  theory  has  been  developed  for 
evaluation  of  the  vibratory  velocity  field  induced  by  an  operating  propeller 
in  both  uniform  and  nonuniform  inflow  fields.  The  analysis  demonstrates  that 
in  the  case  of  nonuniform  flow  the  velocity  at  any  field  point  is  made  up  of 
a  large  number  of  combinations  of  the  frequency  constituents  of  the  loading 
function  with  those  of  the  space  function  (propagation  or  influence  function). 

A  numerical  procedure  has  been  developed  adaptable  to  a  high-speed  digital 
computer  (CDC  6600)  and  the  existing  program,  which  evaluates  the  steady  and 
unsteady  propeller  loadings,  the  resulting  hydrodynamic  forces  and  moments  and 
the  pressure  field,  has  been  extendeo  to  include  evaluation  of  the  velocity 
field  as  well."  TW-*  program  should  thus  become  a  hfghly  versatile  and  useful 
tool  for  the  ship  researcher  or  designer. 
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ABSTRACT 

An  analysis  based  on  the  lifting  surface  theory  has  been  developed 
for  evaluation  of  the  vibratory  velocity  field  induced  by  an  operating 
propeller  In  both  uniform  and  nonuniform  inflow  fields.  The  analysis 
demonstrates  that  In  the  case  of  nonuniform  flow  the  velocity  at  any  field 
point  Is  made  up  of  a  large  number  of  combinations  of  the  frequency 
constituents  of  the  loading  function  with  those  of  the  space  function 
(propagation  or  influence  function).  A  numerical  procedure  has  been  developed 
adaptable  to  a  high-speed  digital  computer  (CDC  6600)  and  the  existing  pro¬ 
gram,  which  evaluates  the  steady  and  unsteady  propeller  loadings,  the 
resulting  hydrodynamic  forces  and  moments  and  the  pressure  field,  has  been 
extended  to  include  evaluation  of  the  velocity  field  as  well.  This  program 
should  thus  become  a  highly  versatile  and  useful  tool  for  the  ship  researcher 
or  designer. 
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NOMENCLATURE 


a  n/u 


f(u) 

9l,2(u) 

h(u) 

K(x) 

m 

K(m) 

j^(m,n) 

x,r,cp 

k 

Up,etf) 

i 

m 

N 

n 

n 

AP(5,P,e0) 

q 

R 

r 

r 

o 

S 


function  defined  in  Eq  (19) 
functions  defined  in  Eq  (18a, b) 
function  defined  in  Eq  (17) 

modified  Bessel  function  of  first  kind  of  order  m  and 
argument  x 

Bessel  function  of  order  n  and  argument  x 
modified  Bessel  function  of  second  kind 
space  function 

derivatives  with  respect  to  x, r,cp,  respectively,  of 
modified  space  function  (after  chordwise  integration) 

variable  of  integration 

blade  loading,  Ib/ft 

integer  multi  pie 

i nteger 

number  of  blades 
blade  index 

order  of  chordwise  modes 

2 

pressure  jump  across  propeller  blade,  Ib/ft 

order  of  harmonic  of  velocity  field 

Descartes  distance 

radial  ordinate  of  field  point 

propeller  radius 

lift! ng  surface 
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t 

U 

u 

*,y,z 


x,  r.cp 


6(x) 


O' 


A(n) (x) 
X 


p 

pf 

a 

T 

$ 

<P 

n 


time,  seconds 

free  stream  velocity 

=  k+a^N  ,  variable  of  integration 

Cartesian  coordinates  of  field  point 

cylindrical  coordinates  of  field  point 

Dirac  delta  function  of  x 

angular  position  of  point  on  propeller  blade,  in  propeller 
plane,  with  respect  to  midchord  line 

angular  chordwise  location  of  point  on  blade 

projected  propeller  semichord  length,  radians 

=  (2n/N)(n-l),n=1,2, ...N 

defined  in  Appendix  A 
order  of  harmonic  of  loading 

cylindrical  coordinates  of  point  on  propeller  blade 

radial  ordinate  of  point  on  propeller  blade 

3 

mass  density  of  fluid,  slugs/ft 
angular  measure  of  skewness 
variable  of  integration 
velocity  potential 
angular  position  of  field  point 
angular  velocity  of  propel ler 
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INTRODUCTION 

The  characteristics  of  the  oscillatory  propel ler-i nduced  velocity 
and  pressure  fields  constitute  the  essentials  for  knowledge  of  the  vibratory 
effects  of  marine  propellers  on  nearby  bodies  (hul 1+appendages) .  Since 
presen'-  technical  developments  put  great  emphasis  on  increased  ship  speed 
these  vibratory  effects  have  become  a  very  Important  problem  to  the  naval 
architect. 

The  basic  tools  in  the  endeavor  to  evaluate  the  pressure  and  velocity 
field  generated  by  the  propel ler  asa  1 ifting  surface  operating  in  nonuniform 
flow  are  found  In  the  fundamental  Davidson  Laboratory  studies^'  ^  where  the 
loading  problem  has  been  treated  by  means  of  the  acceleration  potential 
method.  The  resulting  integral  equation,  which  relates  the  unknown  loading 
distribution  to  the  known  velocity  distribution  normal  to  the  blades,  has 
been  solved  by  using  the  mode  approach  and  collocation  method  in  conjunction 
with  the  so-called  "generalized  lift  operator"  technique. 

The  vibratory  pressure  field  has  been  the  subject  of  a  study^  In 
which  the  analysis  and  numerical  procedure  have  been  developed  and  adapted 
to  a  high  speed  digital  computer.  With  the  knowledge  thus  acquired  of  the 
vibratory  pressure  field  around  an  operating  propeller,  the  vibratory 
effects  on  a  nearby  flat  boundary  can  be  determined  approximately  by  inte¬ 
grating  double  the  free  space  pressure  over  the  flat  surface. 

This  information,  however,  will  not  be  sufficient  for  the  study  of 
vibratory  effects  in  the  case  of  a  body  of  arbitrary  shape,  if  great 
accuracy  is  required.  The  presence  of  a  nearby  body  or  ser’es  of  bodies 
should  be  taken  into  account  by  studying  the  mutual  interaction  between  the 
marine  propeller  and  the  neighboring  bodies.  The  propeller-rudder  interac¬ 
tion  problem  has  been  the  subject  of  a  study^  which  has  considered  the 
complete  interaction  for  any  given  geometry  of  the  interacting  surfaces  when 
both  are  immersed  in  nonuniform  flow.  The  hul 1-propei ler-rudder  interaction 
problem  can  be  tackled  in  similar  fashion  and  thus  the  long  rar>ge  objective 
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of  evaluating  the  propel ler- Induced  vibratory  effects  on  the  hull  is  within 
reach. 

An  essential  requirement  of  such  a  study  Is  detailed  knowledge  of 
the  vibratory  velocity  field.  The  present  study  investigates  the  propeller- 
induced  vibratory  velocity  theoretically  and  on  this  basis  develops  a 
numerical  scheme  adaptable  to  a  nigh  speed  computer  for  evaluating  the 
velocity  field.  Thus  the  existing  program,  which  evaluates  the  steady  and 
unsteady  propeller  loadings,  the  resulting  hydrodynamic  forces  and  moments, 
and  the  pressure  field,  is  extended  to  include  the  velocity  field  as  well. 
This  program  should  become  a  highly  versatile  and  useful  tool  for  the  ship 
researcher  and  ship  designer. 

The  steady  and  unsteady  velocity  field  around  an  operating  propeller, 
in  a  uniform  flow  field  only,  has  been  studied  by  assuming  the  propeller 
either  as  an  actuator  disk^"^,  i.e.  a  propeller  with  infinitely  many 
blades,  or  as  a  lifting  line  model  The  present  study,  an  application 

of  the  unst  ady  propeller  lifting-surface  theory,  takes  into  account  the 
nonuniformity  of  the  inflow  field  as  well  as  the  exact  propeller  geometry 
with  its  finite  number  of  blades  forming  helicoidal  surfaces,  and  hence 
is  unique. 

This  research  was  sponsored  by  the  Office  of  Naval  Research  under 
Contract  N00014-67-A-0202-0022,  NR-062-012. 


VELOCITY  POTENTIAL  DUE  TO  BLADE  LOADING 


The  velocity  potential  at  a  point  (x,r,cp)  due  to  the  loading  of  the 
blades  of  an  N-bladed  propeller  operating  in  a  spatially  varying  inflow 
field  is  given  by 


$(x,r,cf;q) 


WrU  * 


N 

=  11 

n=l  S 


£APW(S,p,e0)e 

\=o 


i\(Qt-e  ) 


<  iXa(r-x)  , 

‘  f  e  ~  i  dxdS 

J  £n  R 
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where  q  and  X  are  positive  integers  defining  the  desired  frequency  of 
the  velocity  field  and  the  appropriate  frequency  of  the  propeller  loading, 
respectively,  and 

x,r,cp  =  cylindrical  coordinates  of  point  In  space  rsrferred 
to  axes  with  origin  at  the  propeller  hub 

§,p,9o  =  cylindrical  coordinates  of  loading  point  on  propeller 

U  =  free  stream  velocity,  ft/sec 

3 

p^.  =  fluid  density,  slugs/ft 

-Cl  =  angular  velocity  of  the  propeller,  rad/sec 
a  -  n/u  ,  1/ft 

AP^(C,p,9  )  =  propeller  loading,  or  pressure  jump  across  the  propeller 
°  biade,  lb/ft2 

R  =  {(t-§)2  +  r2  +  p2  -2rp  cos  j^0o  -  Qt  +  9n  -  a(t-x)-tpj  j-  ,  ft 

“  «  normal  derivative  on  the  helicoidal  surface  at  the  loading  point 
(5,p,9Q)  on  the  propeller 


9n  -  f  (n-0 


n=l,2, . . .N 


S  =  blade  surface,  ft 

The  normal  derivative  to  the  helicoidal  surface  (specified  by  §  =  9^/a) 


(a  ^ _ ^ _ ) 


By  making  use  of  the  expansion  scheme  for  the  reciprocal  of  the 
Descartes  distance  R  ,  vi':.: 


1  =  1 
R  ~ 


I  e 


im(e  -Qt+p  -a(T-x)-q>)  00 


\  I  (lklp)K  (lklr)e 

J  m  m 


(T-§)k 


(when  p  <  ,r  ,  otherwise  p  and  r  are  interchanged  in  the  modified  Bessel 
functions) 
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and  the  facts  that 


N  -l(\-m)0n  rN,  \  -  m  ■  IN  ,  b*  0,  ±1,  ±2,... 
nfj  6  “  M>,  \  -  m  f!  £N 


and  from  the  time-dependent  factor 
\  -  m  «  q  a  j!N 

It  can  be  shown  that  the  potential  function  becomes 

$(x,r,cp;q)  =  JJ  E  AP^(5,p,0  )eIt^tK^m^(x,r,cp;?,p,B  ;q)dS 
S  \=o 

where  the  K-function  is  given  by 


K^(x,r,cp;S,p,eo;q) 


. -Np _ £  e-!mcp 

4n2p,U^ l+a2p2  m«  -«* 
f  m=\-£N 


t  (m)  0 


I •'•“‘'-’h fc  -  hfr  [•  °  I •„< i  - 1 


-oo  p  0  *"  -oo 

It  is  seen  that  theoretically  an  infinite  number  of  frequency  constituents 
of  the  loading  function  AP^  and  of  the  space  function  K  combine  to 
give  the  velocity  potential  at  blade  frequency  and  multiples  thereof,  and 
hence  the  velocity  at  any  field  point.  (This  was  observed  by  Breslin  end 
Tsakonas  in  Ref.  1 1 .) 

The  integration  over  the  blade  surface  can  be  converted  to  integra¬ 
tion  over  its  projection  in  the  propeller  plane,  so  that  Eq  (2)  becomes 


*(x,r,cp;q)  -  j  J  S  APW  (§,p, 9  )e',qntK(m)  — ~-£-  p  dQ  dp 
p  -9b  X=0  °  ap  o 

where  9^  is  the  projected  semichord  at  each  spanwise  location  p  ,  in 
radians.  Since 

=  V  ~  9k  cos9~  = 
o  bcv 

wnere  g  is  angular  position  of  the  midchord  line  of  the  projected  blade 


(4) 


0H 
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from  the  generator  line  through  the  hub,  or  skewness,  and  0  is  angular 
chordwise  location  of  the  loading  point, 


$(x,r,cp;q)  =  J  J  S  L(X)(p,6  )eiq0tK(m)  sin0 id©  dp 

p  o  \=0  aP  a  a 

with  L^(p,0a)  -  APW(§,p,9o).  P0b  ,  lb/ft 

On  taking  the  derivatives  in  (3) 

(m)  e  K(m)  A/l+a2p2  =  +iN  “  Im(0o“9) 


K1 


E  e 

4rrcp,.U  m=-oo 
f  a 


Aia*N('‘-X'f  (.k  -  t)  '(I  k  I  P)K  (|Mr)e'^-5)XdkdT 


The  t- integral  of  (6)  involves 


*  J(ajeN+k)Tj  ,//llllX  -ei(k+a^N)x 
J  e  =  n6(aXN+k) - 


(5) 


(6) 


(7) 


where  6(  )  is  the  Dirac  delta  function.  On  suosututing  Eqs  (7)  and  (4) 
I n  (6) ,  it  becomes 


K'M  =  +  T—Ili  v  im(cr- cp) 

f  m=-co 


{- 


UN(a-ax)  -i*.0bcose< 


a. 


m 


(a  |  £N  |p)  (a  J  £N  !r)[-a2#|+m/p2] 


i  p”  |p)Km(  |k|  r)  (ak+m/p2)  ik(x-a/a)  -i(m-k/a)0  bc°s0 

~  n  J  k+alN  6  e 


1k) 

(8) 


(for  p  <  r  ,  otherwise  p  and  r  are  interchanged  in  the  modified  Bessel 

functions.)  Here  only  pf  and  U  are  dimensional  and  x,p,r,a  and  k  have 

been  nond i mens ional i zed  with  respect  to  propeller  radius  r  as  has  o  in 
/x  o 

Eq  (5). 
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The  loading  function  ^(p>0a)  is  approximated  in  the  chordwise 
direction  by  the  Birnbaum  mode  shapes: 


<M 


*  ”■  {l(X,,)(p)  cot  2 "  +  J  L(X,n)(p)sin(n-l)e  } 

n=2  J 


The  propeller  loading  distribution  at  any  frequency  is  now  given  In  terms  of 

a  spanwise  (p)  loading  distribution  which  is  available  from  the  exlstlnq 
(12  13)  3 

program  *  and  a  chordwise  (e^)  distribution  of  n  modes  which  will 

allow  the  chordwise  integration  to  be  performed  analytically. 

After  the  chordwise  0tf- integration,  Eq  (5)  becomes 


*(x,r,<P;q)  =  J*  J,  S  (p)eiq0tR^m,n)dp 

p  n=l  X=0 


(9) 


where  the  modified  space  function  is 

.if  K'W.e(;)si.«A 

o 

0(1)  »  cot 

®(n)  a  510(0-1)60,  ,  n  >  1 
With  \-m=q«^N^»0,±l,±2  for  each  l 

j^(m,n)  _  INe 


..  im(cr-cp) 

~5r7pfUa~  {e  a  ^  ^  “  a2^N)'m(al^|p)Km(a|i;N|r)A^n^(\eb) 

T  p  w 


-;U>(x-a/a,^^.m(MP)Km(|Mr) 

"00 

*  A(n)((m-|)eb)dk}  (10) 

for  p  <  r  (for  p  >  r,  p  and  r  are  interchanged  in  the  modified  Bessel 
functions).  The  symbol  A^(x)  is  defined  in  Appendix  A. 
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VELOCITY  FIELD 

The  velocities  In  the  x,  y  and  z  directions  of  the  Cartesian  coordi¬ 
nate  system  with  origin  at  the  propeller  hub  at  a  point  x^r,*?  due  to  the 
propeller  operating  In  spatially  varying  inflow  will  be  the  derivatives  of 
the  velocity  potential  $  with  respect  to  x,y  and  z,  respectively.  Since 

z  =  r  cos<p 
y  ®  r  sincp 

the  derivatives  with  respect  to  y  and  z  can  be  obtained  In  terms  of 
derivatives  with  respect  to  the  cylindrical  coordinates  r  and  cp  as 

$  $  $ 

U*  -  sincp  +  -  -Q2  coscp 

z  r  1  w  .  „ 

F  =  f  coscp  "  7  u  smCp 

in  non-dimensional  form. 
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Radial  component: 


~(x,  r,cp;q)  =  J*  _£  j  (p)e1<^tR^m,n^dp 

p  n—  1  \»0  r 


03) 


wh  re  for  p  <  r  for  each  l  and  q 


%p,U  ar 
t  o 


.)  K<"-S>  -  ♦  rialiNL,i  .a^N)e'<»(a-ax),  (.|i|||p) 

'  b*rn  IlS,  ^  ^  111 

‘  CKtn-l(a|iN,r)+'Wl(a,iN  ,r)]  A<n>H) 

+  k  f  1  kfa/0-^  '■<  lkl  P)[Km-l  (I  W  r)] 

*"00 

•  e'k^"^3)^)  ( (rn-k/a)  0b)  dk|  (14a) 


for  p  >r 


b)  ♦  JBS^plg^.^  -AN)eUN^Km(aUN|p). 


4ttp^-U  ar 


[lm.,(a|«  !r>+l)rrt_|  (a  |iN|r)]  hM  (Mj 

•5fM^P\H»)r^,(lb|rHlw,(lb|r)] 

-00  -J 

.e!k(x-a/a)A(F)((m.k/a)%)dk| 


and  for  p  =  r  it  can  be  shown  that 
_(m,n) 

c)  Kr  (p=r)  =  [Eq  14(a)  +  Eq  14(b)]  evaluated  at  p=r 


Tangential  component: 

*®(x*ri9;q)  ,  „  ®  •  /,  -x  .rs.  ,  -x 

-J— - -  =  1  f  J  V  L^n^o'»e,q^tK^ni,n^- 

r  r^U  r  Jp  ^0L  (P>e  Kcp  * 


04b) 


(14c) 


(15) 


where  r  is  nondimens ional  with  respect  to  rQ  propeller  radius  and, 
for  p  sr  and  each  jJ  or  q 
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gOn,n) 


+ 


o 

4ttp-U  ar 
f  o 


l„(a|«N|p)Ktn(a|itN|r)A(;')(X9b) 


_  j.  p”(ak+m/p^) 

IT  J  k+aiN 
-00 


'm(|Mp)K|n(|k|r)clk(x-'r/a) 


•  A^Km-k/aJg^dkj-  (16) 

For  p  >  r,  p  and  r  are  Interchanged  In  the  modified  Bessel  functions. 

The  Integral  terms  of  the  modified  kernels  (Eqs  12,  14,  16)  all  have 
Cauchy-type  singularities  at  variable  points  depending  on  values  of 

<=°’:tl’±S .  Use  of  a  transformation  u-k+afN  fixes  the  singular  point 

at  the  origin,  thus  avoiding  the  complications  arising  In  the  numerical 
integration  from  singular  points  at  varying  values  of  k 


influence  Functions 

Letting  u-k+a^N  ,  Eq  (12)  becomes 

i?(m,n)_  Ne’ta-'n'P-aWx) 

x  ”  !  2 

%p  U  ar 
r  o 

•  {a«N(a2«N-m/p2)lm(a|fN|p)Km(aUN  |r)A(">(Xeb) 

.  1  f  h(u)-h(-u)  ,  1  , 

n  J  u  'duf  for  p  ^  r 
o  J 

where  h (u) « (u-a^N)  (au-a^N+m/p2)  lm((u-aXN  |P)Km(|  u-a^Nj  r } e  1  u (x-cr/a) 
•  A(n)((\-u/a)eb) 

The  integrand  at  u=0  is  evaluated  in  Appendix  B  as 
a)  when  q=i=0  ,  m=\=0 


07) 


,ifnrMHl.-h(-u)1 
U-aO  L  U  J 
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b)  When  q«A»0, 

]  -  r(l/p2)(p/r)mA(n)(\0b)  for  p  ^  r 

l(Vp2)(r/p)mA(n)(\0b)  for  par 

c)  When  q*ANi<0 


1 1m 
u-*0 


+ 


.  ,2  |i„|p)Km(a|  ANtr)} 

•  [[  2a2AN-m/ p2- 1  (x-<j/ a)  a-^N  (a2AN-m/p2)  ]A^  (A©b) 
-l0b£N(a2^m/p2)A|n)  (\0b)] 

{pKm(a|  AN  |  r)[lm_,(a  |AN  [p+l^,  (a  Un  |  p)  ]-rIm(a  I  AN  ;p) 

’  [V]^al1Nlr^+KhH-l^aliNl  r)]}‘  ■JW(a2AN-m/o2)A^n^(\0b)  ,  for  p*r 


Where  the  upper  sign  (-)  Is  used  when  A  >  0  and  the  lower  sign  (+)  when 
A  <  0  .  When  p  >  r  ,  p  and  r  are  interchanged  ir  the  factors  In  braces 
(Involving  modified  Bessel  functions  and  their  derivatives). 

The  closed  term  of  (17)  can  be  shown  to  be  equal  to  zero  when  q*A«=0, 

Equation  (14a)  for  p  <  r  becomes  with  u-k+aAN 


i  .  .2 


p 

4rrprl)  ar 
f  o 


•  {-  |AN  |(a2AN-m/p2)lm(a|AN|p)[Krti_1(a|AN|r)+Km+I(a|AN|r)]  A(n)(\0b) 
j  “  g,(u)-g,(-u) 


+  2tt  I 


(l  8a; 


wh«re 


9 ,  (u)  =|  u-aAN |  (au-a2 i-N+m/p2)  lm(|  u-aAN  |p) [k^_ ^  ( |u-a  N  |r)+K^,(|  u-aAN  |r)] 

.eiu(*-a/a)A(S)((Wa)eb) 

It  can  be  easily  shown  that  when  q=A=0  the  closed  term  of  (18a)  becomes 
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(1/2)  (p/r)m(m/rp2)A^  (\6^) 
where  m**\  since  q=0. 

The  Integrand  at  u-0  is  evaluated  in  Appendix  C  as 

a)  When  q*£*0,  m=\=0 

rg.(u)-g,(-u)-, 

,  }m  [ J - - - 1  =(4Vr)Aln)  (0) 

u  -0L  u  J 

b)  When  q=i»0,  m=X/0 

rg.  (u)-g,  (-u)-,  ,  ,-x 

1 } m  - - - J  =  (2/-) (p/r) mj[a-t{im/f£) (x-cr/ a)  ]a  (X0b) 

+  ( im/p2)  (eb/a)A^  (A6b)| 

c)  When  q=/N/0 

lim  [il— ii_]=  ±sfl  (a|«|p)rK  ^.IWIIrJ+^^liNlrT} 

u  -0  L  ^ 

•  |J^2a2AN-m/p2-i  (x-o/a)  (a£N)  (a2^N-m/p2)J  A^  (\0b) 

-?0b/N(a2XN-m/32)A{n^  (X9b)} 

+aAN(a2A  N-m/p2)A^ (X0b) 

*  |  r)+K(|H_1  (a  |  jBN|  r)J 

-r  lm(a  | AN |  p)[Km_2(a  \M  j  r)+2Kn(a |  iH  |  rJ+K^a  |£N  | r)  J- 

where  the  upper  sign  (+)  is  used  when  l  >  0  and  the  lower  sign  (-) 
when  l  <  0  , 

Equation  (14b)  for  p  >  *  becomes  with  u=k+a^N 


1 1 
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-fm  n)  M  i  (X<7-tlKp-aiNx) 

It)  ’"'(p  >  r)  -  -  -N-S-  - 

%p£lr ar 
f  o 

.  |-4<l/2)a  |^H|(a2AN-m/p2)Km(aj£N|p)^lm  l(a  |iN  (rj  +  t^ j (a|  jtN|  r)  Ja^  (X6b) 


(l/2rr)  J 


“  g2(u)-g2(-u) 


-  du} 


08b) 


-ce 


“here  92(u)=  -  lu-aWllfau-a^N+m/p^K^lu-aiN  |p)[lm]  (|u-a«|  rj+l^,  (]u-atN|  r)] 

•  eiu(''-e/a)A<")((x.u/a)0b) 

In  this  case  when  q=£=0  the  closed  term  is 

-(i/2) ( r/p) m (m/ rp2) A  ^ ^ (X0b) 

The  integrand  at  u=0  is  evaluated  in  Appendix  C  as 
a)  when  q=/=0  ,  m=X*=0 


1  im 

u  -*0 


[■ 


g.j(u)-g0(-u) 


u 


b)  when  q=£=0  ,  m=X/0 
-g2(u)-g2(-u) 


I  im 
u  -*C 


P 


c)  when  q=AN^0 

-92(u)-g2(-u) 


1  im 
u  -0 


P 


■J  =  (-2/r)(r/p)m|[a+(im/p^(x-cr/a)]A^  (X9b) 

+  ( im/p2)  (eb/a)A|n^  (X9b)} 

■]  =  +2Km(a|  ^lP[lrn.1(a|iN|r)+lmfl(a|£N  |r)] 

•  •j2a2JJN-m/p2-i(x-a/^)ai5N(a2/N-m/p2)jA^n^  (x©b) 

-ieb/N(a2/N-m/p2)Ajn' (X0b)} 
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+(a£N)  (a2£N-m/p2)A^n-  (X6b) 

*  {p[Km.l<aiAM|p)+Km+l<al^lp)‘I,m-l<a  1“  lr)+,«H-,  (*l  "I  r)] 

-rKm(a  iiN|p)[lm-2(a!i.N|r)+2lm(ai£N|r)+»^2(a|£N|r)]} 

where  the  upper  sign  (-)  is  used  when  £  >  0  and  the  lower  sign  (+) 
when  £  <  0  . 

With  u*k+a£N  ,  Eq  (14c)  for  p«r  becomes 
K^m,n^ (p*r)  [Eq  (18a)  +  Eq  (I8b)j  evaluated  at  p*r  (18c) 

On  substituting  u=«k+a£N  Eq  (16)  becomes 

j;(tn,R)  _  _  ^Ifro-nT-alllx) 

^  4np  u^ar 

f  o 

•  {m(a2£N-m/p2)lm(al£N  |p)Km(a |£N  | r)A^  (\0b) 


where  f  (a)  - (au-a2£N4m/p2)  i m ( |  u-a£N  | p)  Km ( |  u-a£N  |  r)  e 1 U (  (\-u/a)  0b) 


it  can  be  shown  that  when  q«£=P  the  closed  term  of  (19)  becomes 
(-m/2p2)(p/r)mA^  (X.0^)  for  p  s:  r 

(r/p)mA^ (X9b)  tor  p  ^  r 


The  integrand  of  (19)  is  evaluated  for  u*°0  in  Appendix  D  as 
a)  when  q“£»0  ,  m=\=0 


1 1m  m 
u  -*0 


0 
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b)  when  q*A*=0  ,  m=Xi*0 

lim  ra  piiillfllldj  s  jj^a+j^-e/^m/p2^^  (X0b)+(i0b/aj  (ri/p2)Aj^  (X0b)| 

•  f  (p/r)m  for  p  £  r  ^ 

\(r/p)m  for  p  ^  r  J 

c)  when  q=/N^0  ,  m=X-AN 

1  im  m  [iMl£LuI]  =:  +aJT a-i  (x-<j/a)  (a2AN-m/p“)"|Ain^  (X0b) 
u  _.0  L  u  -*  m_  J 

-(?0b£)(a2iN-m/p2)A|n)  <^©b)}*m(a  [AN  jpKja  |AN  |r) 

-m(a2AN-m/p2)A^  (X0b) 

•  |pKm(aUN|  r)J^lfn_1  (aj  AN  |p)  +  lfl1+1  (a |  AN  |p)J 
-rlm(a|AN|p)[Km_1(ajAN|r)+Km+1(a  |AN|r)]j- 

for  p  £  r  ,  otherwise  p  and  r  are  interchanged  in  the  product  of  the 
modified  Bessel  functions  and  its  derivative  (the  factor  in  braces  in  the 
second  term).  Here  the  upper  sign  (+)  is  used  when  A  >  0  and  the  lower 
s  i  gn  (-)  when  A  <  0 

When  m=0,  i.e.  X=AN  whate  'er  the  value  of  A  , 

R(m,")  -  0 


STEADY-STATE  AND  BLADE-FREQUENCY  VELOCITY  FIELD 
FOR  NONUNIFORM  AND  UNIFORM  INFLOWS 

It  is  easily  seen  from  Eq  (9)  that  in  the  steady-state,  nonuniform 
inflow  case,  when  q=A=0  ,  the  velocity  potential  is  given  by 


$(x,r,<P;0)  »  R.P.J  _E  E  L(?v,n)  (p)  K(^n)dp 
p  n= 1  X=0 


(20) 
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:  l 


i  1 

t 

i 

]  l 
!  I 
I 
I 


whereas  in  the  unsteady  blade-frequency  case  where  |jjj*  l  the  velocity 
potential  becomes 

i(x.r,<p;N)  *  R.P.  {e’^J  S  s{L^,n^(p)  K^'N,n) 

L  p  n=l  X=0C 


conj  .  [|>>5>(P)  K^’Jdp} 


where  in  the  first  term  on  the  R.H.,  i5  evaluated  at  m=*X-N,A*l  , 

and  in  the  second  term  at  m=\+H,i=-l  .  Therefore  the  blade-frequency 
velocity  potential  is 

*(x,r,o);N)  -  j  £  Z  (JX,n*(p)  R^'^+conj  (p)R*X+N,n*"i}dp 

p  n=l  \=0  1  L  M 

(21) 

In  t'  uniform  inflow  case,  on  the  other  hand,  with  \=0,  the  corres¬ 
ponding  values  are: 


and 


$(x,r,cp;0) »  J*  _£  L^0,n^(p)  K*0,n^dp 
p  n=I 

$(x,r,cp;N)=  J  £  (p)["R  +  conj  .R(N,n^"]c!p 

p  n=l  L  J 


(22) 


r 


lr 

i 


=  2  J  _£  L^°,n^(p)  R^'^dp 

p  n=l 


(23) 


The  velocity  field  is  given  by  the  derivatives  with  respect  to  x,y,z 
of  these  velocity  potentials.  In  non-dimensional  form  in  the  steady-state, 
nonuniform  inflow  case,  q=5=0  ,  m=\ 

v!0)  *„<»> 


v(0>  $  (o) 


u  u 


-  J  J  E 

l>’">(p)  K(X’S)dp 

X 

(24a) 

p  n=l  X=0 

=  J  S 

L<X,n)  (p)  j^(x^)dp 

(24b) 

p  n=l  X=0 
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V<0)  #„(0) 


-f-  -  tVu  ■  ?I  -s,  E  L(X'n)(p)  K<,X>n)dp 

o  p  n**1  \=0  T 


(24c) 


where  j^(m>n)  an(j  are  g|ven  by  £qS  (J7),  (18a,  b  or  c) 

and  (19),  respectively.  Note  that  the  argument  of  the  potential  refers  to 
the  value  of  the  frequency. 

In  the  unsteady,  blade-frequency  case  q=N,  |x|  =1,  m=X-AN 
Vf  #„(N) 


u 


=  — - =J  _£  E  •fL^n^(p)K^"N,n^ 

p  n=l  \=0  ^  A 


+conj  (p)R^+t,,n)]}dp  (25a) 


v‘N)  i  (N) 


“  =  “  '  ■''p  J,  Jo  (L  X'n 


«onj.[L^S>(p)K<W“’S)]}dp 
*m(N)  .  r  t\  “\  u(^-N,n) 

^‘^^-E,JoL  <p),p 

"  f-  n=l  X=0 


,(N) 


(25b) 


+conj  .[l  n)  (p)  K^+N' n)  ]}dp  (25c) 


NUMERICAL  RESULTS 

A  numerical  procedure  for  the  evaluation  of  the  velocity  field 
induced  by  a  propeller  operating  in  spatially  uniform  or  nonuniform  inflow 
has  been  developed  and  incorporated  in  the  existing  program  which  computes 
the  loading  on  the  interacting  blades,  the  resulting  steady  and  unsteady 
hydrodynamic  forces  and  moments  and  the  pressure  field. 

The  enlarged  program  has  been  used  to  calculate  the  steady-state  and 
blade-frequency  theoretical  velocities  in  the  x,  y,  and  z  directions  in  the 
neighborhood  of  the  NSROC  3-blade,  12-inch  diameter,  marine  propeller 
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No.  4110  of  blade  area  ratio  BAR  **  0.6,  which  had  been  treated  In  earlier 
studies.  The  blade  loading  distributions  and  resulting  forces  and  moments 
had  been  presented  in  Refs  3,  5  and  12,  the  blade-frequency  pressures  in 
the  vicinity  of  the  propeller  in  Ref  6. 

In  Ref  6  the  theoretical  pressure  calculations  were  shown  to  be  In 

(14) 

satisfactory  agreement  with  the  available  experimental  results'  '  for  both 
uniform  flow  and  the  nonuni  form  flow  generated  by  a  wake  screen.  The 
present  computations  of  the  velocity  field  have  been  performed  for  the  same 
conditions,  viz., 

Advance  ratio,  J  0.833  (Design) 

'ip  Clearance  5%  of  diameter  (r  =  1.1) 

Fngle  cp  (clockwise  from  0 

i2  M  position  looking  aft) 

over  a  range  of  ax>al  distance  x  from  the  propeller  plane,  from  -0.3  of 
radius  (upstream)  to  -+0.3  radius  (downstream).  The  harmonic  analysis  of  the 
screen  wake  survey  was  that  supplied  by  NSRDC. 

The  results  of  the  calculations  are  shown  in  the  figures.  Figures  1 
and  2  present  the  real  and  imaginary  parts  of  the  unsteady,  blade-frequency 
x,  y  and  z  components  of  the  velocity,  in  nonuniform  and  uniform  flow, 
respectively.  Figures  3  and  4  depict  the  steady-state  components. 

Figures  5~7  are  comparisons  of  the  unsteady,  blade-frequency  veloci¬ 
ties  for  uniform  and  nonuniform  flow  conditions.  These  figures  show  that 
nonuniformity  in  the  flow  reduces  the  amplitudes  of  the  velocities  in  the 
neighborhood  of  the  propeller.  This  was  to  be  expected  from  the  comparisons 
of  the  experimental  blade-rate  pressures  for  the  same  tip  clearance,  angle  cp 
and  J  in  nonuniform  and  uniform  flow  given  in  rigs  33-35  of  Ref  14,  where 
it  is  seen  that  nonuniformity  in  the  flow  increases  the  pressure  magnitudes.* 

The  pressure  can  be  computed  from  the  velocity  components  by  using 
Bernoulli's  linearized  equation  for  free-space  pressure 


*fc 

It  is  to  be  noted  that  the  peak  in  olade-rate  pressure  in  uniform  flow  in 
Fig  34  of  Ref  14  is  in  error,  as  comparison  with  Figs  15  and  16  of  that 
reference  wi 1 1  show. 
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PtK,r,<p;t)  -  -Of<||  +  U  |J) 

■  -pf(flrr  V_  +  UV  ) 
f'  09  x 

(The  pressure  on  a  flat  plate  Is  twice  the  free-space  value.)  Calculations 
of  the  blade-frequency  pressure  In  a  uniform  flow  field  for  the  small  axial 
range  -0.3  <  x  <  0.3,  using  the  velocity  components  shewn  In  Figure  2,  check 
quite  well  with  the  corresponding  pressures  on  a  flat  plate  arising  from 
blade  loading  in  Figure  13  of  Reference  6,  which  compared  well  with  experi¬ 
mentally  obtained  values. 

CONCLUSION 

The  vibratory  velocity  field  induced  by  an  operating  propeller  In 
uniform  and  nonuniform  inflow  fields  has  been  determined  on  the  basis  of 
lifting  surface  theory,  and  the  developed  numerical  procedure  has  been 
incorporated  in  the  existing  program,  adapted  to  the  CDC  6600  high-speed 
digital  computer,  which  evaluates  the  steady  and  unsteady  propeller  load¬ 
ings,  the  resulting  hydrodynamic  forces  and  moments,  and  the  pressure  field. 

The  analysis  demonstrates  that,  in  contrast  to  the  case  of  uniform 
flow,  in  nonuniform  flow  the  velocity  at  any  field  point  is  made  up  of  a 
large  number  of  combinations  of  the  frequency  constituents  of  the  loading 
function  with  those  of  the  space  function. 

The  results  of  calculations  for  a  3-blade  propeller  in  uniform  flow 
and  in  a  3-cycle  screen-generated  v  „ke  show  that  the  effect  of  the  non¬ 
uniformity  in  the  flow  is  to  reduce  the  amplitudes  of  the  velocities  in  the 
neighborhood  of  the  propeller.  Although  no  experimental  velocity  data  are 
available  for  this  case,  this  effect  of  the  nonuni formi ty  is  confirmed  by 
the  experimental  data  on  pressure  fields  induced  by  the  same  propeller 
under  the  same  conditions  which  show  that  the  nonuniformity  increases  the 
pressure  amplitudes.  Furthermore,  pressures  computed  from  the  theoretical 
velocity  components  in  a  uniform  flow  field  show  good  agreement  with  both 
earlier  ca'culations  of  the  pressures  due  to  blade  loading  and  experimental 
pressure  measurements. 
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The  particular  case  of  a  propeller  in  a  screen  wake  was  selectee  for 
purposes  of  comparison  and  because,  the  wake  being  pure  third  harmonic  and 
strong,  wake  measurements  were  accurate.  Unfortunately  these  wake  charac¬ 
teristics  are  not  those  of  a  typical  ship. 

It  Is  to  be  noted,  furthermore,  that  the  present  report  has  con¬ 
sidered  orly  the  effect  of  blade  loading  on  the  velocity  field.  Blade 
thickness  Is  also  important.  An  expression  for  the  velocity  field  due  to 
this  effect  has  been  developed  analytically,  and  will  be  Incorporated  In 
the  present  p.ogram  and  applied  to  the  case  of  a  propeller  In  an  actual 
ship  wake,  In  an  extension  of  this  Investigation. 
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AXIAL  DISTANCE  X,  RADII 


BLADE-FREQUENCY  VELOCITY  FIELD  AROUND  PROPELLER  4118  IN 
NONUNIFORM  INFLOW  (  3-CYCLE  SCREEN  WAKE) 
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AXIAL  DISTANCE  X,  RADII 


FIG, 2.  SLADE-FREOJENCY  VELOCITY  FIELD  AROUND  PROPELLER  4118  IN 
UNIFORM  INFLOW 
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FIG. 3.  STEADY-STATE  VELOCITY  FIELD  AROUND  PROPELLER  4118  IN  NON- 
UNIFORM  INFLOW  (3-CYCLE  SCREEN  WAKE) 
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FIG.  4,  STEADY-STATE  VELOCITY  FIELD  AROUND  PROPELLER  4118  IN 
UNIFORM  FLOW 
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F!G. 5.  UNSTEADY  BLADE-FREQUENCY  Vx /U  VERSUS  AXIAL  DISTANCE 
X/r0  FOR  UNIFORM  AND  NONUNIFORM  FLOW  CASES 
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FIG. 6.  UNSTEADY  BLADE  FREQUENCY  Vy/U  VERSUS  AXIAL  DISTANCE 
X/r0  FOR  UNIFORM  AND  NONUNIFORM  FLOW  CASES 


FIG, 7.  UNSTEADY  BLADE-FREQUENCY  Vz/U  VERSUS  AXIAL  DISTANCE 
X/r0  FOR  UNIFORM  AND  NONUNIFORM  FLOW  CASES 
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APPENDIX  A 
A^(x)  and  aJ^  (x) 

•  A^(x)  ”  “  J  0(n)  expt-lxcosQ^Jsine^  d0a  (A-1) 

o 

where  9(1)  =  cot  (0^/2) 

9(11)  »  s  i  n  (n- 1 )  0^  ,  n  >  1 

Then 

A(0(x)  -  Jo(x)  -IJ,(x) 

—  n-2 

A(n)(x)  .  Izll -  [j-_2(x)  +  J-(x)J  ,  n  >  1 

where  Jn(x)  is  the  Bessel  function  of  the  first  kind  of  order  n  and 
argument  x 

II  aJ^(x)  J  ®(n)  exp (  - 1 xcos 0^) cos e^s  i  0^  (A-2) 

o 

a}°(x>  •  \  [J0<x)  ■  J2^x0  *'J|M 

Ain>  (x)  "  [Jn-3(x)  '  Jn+l(x)]  '  "  >  1 


It  Is  to  be  noted  that 


(n)  ..  j  tt 

e  “  ^  J  0(n)  exp(-ixccsea)coseasine^ 


de 


a 


-  -Ia{5)(x) 


A-1 
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APPENDIX  8 

Evaluation  of  the  Integrand  of  at  u  ■  0 
The  integrand  of  Js 

j-h(u)  -  h(-u)1 

L 

where  h(u)  -  (u  -  aiN)  (au  -  ?2£N  +  -^)  ^ln|(lu-ajtNlp)Km(lu-ajjNlr)^ 

•  1>  a<">((x-  f)  eb) 


(B-l) 


and  q  ■  !H,  l  ■  0,  ±1,  ±2,... 

I.  q  -  l  ■  0,  m  «  > 

h(u)  =  u(au+  -j)  rim(up)Km(ur)"J  e,U*x_  a*  A^((\-  £)eb) 

P 

h(-u)  -  -u(-au+-^)  [lm(up)Km(ur)  j  e",u^x“  a*  A^((\+  £)eb) 

p 

Then 


■  jJm(up)Km(ur)]  fau[e,,(x-|)  A(n)(U- j)eb) 

.e-ia(x-|)  A<"»C(x+  ^e„)T 

+  4[elu'x-  f>  A(">((X-  f)eb)^_‘U<X-  A(S,«X+  >b>l} 


a)  when  m  *  \  t  0 


m 


lim  [Mu)-h(.-u)~|  JP 


j  (f)  (X0b)  for  p  r 


u-*0 


\  (“)  A^(xeb)  for  p  >  r 

P  P 


(B-la) 


B-l 


R-1588 


b)  when  m  ■  \  «  0 


11m  n 

u-0  L 


h(u)  -  h(-u)i 
u 


(B— l  b) 


II.  q  »  £N> 0 


Hm  .  1  lm  i  |(u-ajtN)  (au-aZ£N+  ~2)  [lm((ajtN-u)p)Km((ajtN-u)r)  j 

P 

•  elu(x-  f)  A<S)((>.-  f)0b) 

-  (-u-a/N)  (-sm  -aZm~)  jjm((ajtN4-u)p)Km((a£N+u)r)J 


a  A  0# ( (X+  ^)eb)}  for  p  <  r 


By  L'HospItal's  rule 


I  im  rMu).-h(-u?1  -  SlilMiMiiil] 

u-0  L  u  J  du  VO 


1 1m  Ij^au-a^ihH-  -y-  +  au-a2£N+l(x-  |)  (u-aXN)  (au-a2|N+  -2)~j 
u-»0  p  p 

•  [l„((a«l->i)p)Km((alN.U)r)]  ,lu(x-  |»A(S)((  } 

3  b 

+  |^-au-a2xWf  *y  -  au  -  a2£N  +  f  '  £) (-u-aiN) (-au-a2£N+  *y)l 

p*  p 

•  [«m((ajeN+u)p)Km((axN+u)r)J  e_iu(x"  aV((\+  £)eb) 

+  (u-ajjN)  (au-a2£N+  -y)eiu^x“  ((\-  i-)tf-plrJ,((aJeN-u)p)Km((aJeN-u)r) 

P  L 

-  r  lm((a£N-u)p)Kj)((a,eN-u)r)] 

-  (-u-ajeN)(-au-a2£^  2L2)e”,u(x-  f)A(n)((\+  £)eb)  [pi;((a^u)p)Km((a^u)r) 


+  r  lm((a^W+u)p)K^  ( (aj^N+u) r) 


B-2 
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n  m  r  -1  lufx-  2*'  ^((X“"a)9i.) 

+  (u-aiN)(au-a2xN+ -j)  [lm((ajtN-u)p)Km((ajeN-u)r)j  e  a'  -  1 — - 

7  m  r  -1  .|„/u  2A<^  *((X+ 

-  (-u-ajtN)  (-au-aZjeWf  -^)  [lm((aiN+u)p)Km((a^N+u)r)J  e  lu'x"  a> - - - •’ 


where 


.tl* 

3u  I  ,._n  a 


*l^(X9b) 


a«  *u=o  a  ‘  b 

and  (XS^)  Is  defined  In  Appendix  A. 
Then,  whatever  m, 


Hr,  [■ 

U-»0  L 


^ . "  “  2  [lm(aUNIp)  Km(aljgNlr)] 

.  {[+2a2iN-  -y  -  I  (x-  J)  (a£N)  (a2^-  ^)]  A(n)  (X0b) 
P  P 

-  i  ebjeN(a2XN — ^)A{n^  (\0b)]“ 

p 

-  2  (a£N)  (a24N-  ~-)A(n) 

P 

.  [pl(|}(aUNIp)Km(aUNIr)+rlm(aUNIp)KtJi(aUNIr)] 

(B-2) 


where 


■i[|m.i<alw,p>+|Wi(a'1NI»)] 

K^aljJNIr)  -  -  i  [Km_,  (aUNI  r)+Km+j  (al  £NI  rj 


(B-?a) 


Here  p  <  r  otherwise  p  and  r  are  interchanged  in  the  modified  Bessel 
functions  and  the  last  term  in  brackets  of  (8-2). 


8-3 


Ol|C 
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Ml.  q  *  jfcN  <  0,  whatever  m, 


(iJi.-.h.(--.u)r|  ,  Mm  1  {(u-aiN)(au-a2/H4-  -^)  [l  J(u+al  jtNOp^UiH-al  XNl)  rl 

*  *'"<*-$  A(S)((x-^)eb) 

-  (-u-aXN)  (~au-aZjiHf -^)  jjm( (-in-al  jyj|)p)Km( (-u+aUNl)  ' )] 

P 

.  .-'“<*-5)  a(5)<u+X>} 

.  2  [^(.IW'pjK^a'W'r)] 

.  {[2a2jiN-  Jy  -  I  (x-  £)  (a2N) (a2«l-  -y(]  ^  (\\) 

P  P 

-  lebJtN(aZlN-  ^A, (n)  (xeb)} 

p 

+  2(aiN)(aZjlN-  JL)k^  (X\) 

P 

.  [pM(alANIp)Km(al£Nlr)+rlm(a|JeNlp)K||)(aljeNlr)]  (B-3) 


for  p  <  r.  Again  If  p  >  r,  p  and  r  are  Interchanged  In  the  modified  Bessel 
functions  and  the  last  term  In  brackets  in  (B-3). 
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PEND IX  C 


Evaluation  of  the  Integrand  of  R^  at  u»0 
A.  For  the  case:  p  <  r 

The  Integrand  of  R^m,n^(p  <  r) 
r9|(»)-9I(-u)^ 

L  u  J 


(C-l) 


where  g;(u)  - |u-a£N|  (au-a2iN+m/p2)|lm( |u-a^N 

'  ,  ( I u-a^N  |r) H-K^. ,  ( | u-a^H|  r)}e ,u<x-»/a)A(")  ((Wa)eb) 

and  q-/N  ,  . . . 


q-jt-0  ,  m»\ 


9|(u).U(.wp2){,mM[V|(ur)tV|(uo]}a'“<^)/")((x.u/a,v 

a"d  9'  <ur>+KmH(-)]}e-'u<>'-^)A(")((Wa)  %) 

(B-2) 

Then 


[' 


9j (u)-g, (-u) 
u 


<Ur)+V !  <ur)  }{a[e'  u(x-a/aV S)  < fc-u/a)  9b) 

+  e-lu(x-,/a)A(S)((Wa)6b)-j 

*Tp“M»(:|l(Wl)V 

ii  D 


P  u 


1  u  (x/s/a)  A  (")  ( (,i+u/a) 


R-J588 


a)  when  nr=xyo 

9i(u)_g.(“u)  r  _jii  a  r  t~S 

Hm  -J - I -  a  f£|  1  /oA(n), 

u  -*  0  u 


f]  fKW 


*  i< h/p*),ta  ((Mu/a)8b) 

J  u  -0  L  u 


But 


, ,m  |-e"M>fl(IW/.) 9b) Iu ' <x-°/a)/, ,S> ( (X+u/a) e„) 
mL  u  ]- 


u  -*0 


IW.)[2  A<S>(X90)]  ♦  .  ?A(n)((Wa)9b)j  ^ 


since 


with  a|  ^ (x)  as  defined  in  Appendix  A 
Therefore  the  Integrand  at  u=0  is 

7  [?T{[a+^m/p^x"ff/a)]A<n)(xeb)+(im/p2)(0b/a>Ain)(xeb)} 

b)  when  m«=X=0  the  Integrand  at  u=0  is  easily  seen  to  be 

~A(n)(0) 

I  I  q— j0N  >  0 

-9,(u)-g}(-u) 


2l(x-CT/a)A^n^(\eb)  +  21  (eb/a)Afn^(X9b) 

^A^n\(X-u/a)  eb)  I  ( (X+u/a)  0  )  I  - 

9u  u»0=  ?iu  nan*s^®h^)Ai  (**0 


(C—  1  a) 


(C-lb) 


1  im 
u-*0 


1-3)  Vu/  3}  \  «/  -1  |  ,  p 

L - u - J=  limou{(a^N-u)(au-a  XN-Hn/p2)|lm ( (aiN-u)p) 

•  [Km_ ,  ( (aAN-u)  ,  ( (a*N-u)  r)  ]}e 1  u  (x'a/a) A  (n)  ( (x_u/a)  ^ 


C-2 
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-  (a/N+u)  {-au-a2jfcN4m/p2)  -^1  m  ( (aiN+u)  p) 

•  [K,_ ,  ( (a/W+u)  rJ+K^ ,  ( (aAN+u) r)]}e" ' U ^/a\^  ( (X+u/a) eb)} 

0 

0 


By  L'Hospita) 's  rule 

r9.(u)-gi(-u)  paCg.t^-g^-u)! 

u'^oL  u  J  L  Ju*0 

Then,  whatever  m  , 

J  m  (a  |  iN  |  p)^. ,  (a  |  AN  j  r )  ^ ,  (a  |  i  N  |  r)  ]} 

•  -|[2a2jfcN-m/p2-l  (x-cr/a)  (a^N)  (a2jIN-m/p2)  ]A^  (A0b) 
-«ebJeN(a2iN-m/p2)Aj^(\e  b)j- 

+2(aAN) (a2j&N-m/p2)A^ (\0b) 


•  {pP  (a  |  JIN  |p)[Km_,  (a|  JIN  lO+K^,  (a|  |r)  ] 


+  r  1  m (a)  JIN  |p)[K^_j  (a|  £N|  O+K^,  (a|  4N|  r) 
where  r(a|j!N|p)  =  |  (a| -CN  |p  +  1^,  (a  |  AN  |p)  ] 

k;_i  (a  \m\ 0+K^,  (a  |J5N|  r)  -  ~  [K^a  |jtN  |  r)+2Km(a[  IH\ O+K^al £N|  r)] 


(C-2) 


(C-2a) 


III  q=4N  <  0  ,  whatever  rn  , 


r9 1 (u) -g ! (-u) 

L - 3 - ] 


u  -•  0 


=  1  im  ~-j  (u-a^N)  (au-a2^N-Hn/p‘'){ lm((u+a  |j£N  |)p)j" K^_1  ( (u+a  |  ^.N  j) r) +Km+  ^  ((u+aj  £N|)r)l| 
u  -.0  1  ^ 


,iu(x-a/a)A(-n)((N.u/a)eb) 


'>3 
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-(-U-.MO (-.u-a2iH«./p2){l„((-^«<l  )P)[V.((-U+SI‘N1  °]} 

e- Ju(x-a/a)  A(n)  ( (x+u/a)  eb)  j. 

.  -2{lm(a  |M Ip)^, (a|tH  |r)+K^,(»|»  |r)| 

.  |[2a2<N-m/p2- 1  (x-p/a)  (alN)  (a2iN-m/p2)  lA*"'  (>-«b) 

-I  @bWi(a2j(N-fli/p2)A|^  (A0fe)} 

+2(,^)(aWp2)A(5MMb){p';(a|^lp)[V.<alWlr)^l(aUN|r)] 

♦rlm(.|iN  |p)[K;.,(aUN|  r)+K^,  (a  |iN  |r)]f  (C-3) 

B.  For  the  cas*:  P  >  r 

The  Integrand  of  K^.m,^(p  >  r)  is 

^g2(u)“9g(-u) j  (C-4) 

where  g2(u)  is  given  in  Eq  (18b) 

I  q«Jt*0  ,  m*=X 

[-9"(U)u-(-~]a  -uKm(up)[,m-l(ur)+,nH-l(lir)J 

.  u  (x"°'/a)  ( (\-u/ a)  9b)+e"* '  °  (x-or/a  ?  A(  '^  ( (X+u/a)  8fa)  ^ 

+  jl.  Mu(x-a/a)Atf)(^^ 

p2u  L 

a)  when  m=X»<0  it  can  be  shown  that 

rfgj^sLl'L  .(r/p)"(3/r){^(lm/p2)(x-p/a)3A<r,)fr8b) 
u  -»0*-  U 

f  ( ;m/p2)  (8^/9)^! 

C-4 


(C-4a) 
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b)  when  m*\-0 

rgP(u)-go(-u)-, 

llm  f— - - - ■]  -  o 

u  -K)  *■  u 

1 1  q-iN  >  0 

J!m  r-92(u)“92("u)-|  ,,  J.  faXN-u)(au-a2jtN+m/p2)K  ((aih-ujp' 
u-oL  u  -I  *j-<0  u 

•[lm_1((a£N-u)r)+l|^1((aJeN-u)r)]e,u(x“CT/a)A(n)((\-u/a)0b) 
+(alM+u)  (~au-a2AffHn/p2)  Kffl(  (aiN+u)p) 

-  -2Km(SXNp)[lm_J  (aJeNr)+1m+1  (aiNr)] 

|[2a2lN-m/p2-l (x-a/a)aAN(a2AN-m/p2)]A(n) (\0b) 
-!0b4N(a2*N-m/p2)Ajn) (X0b)} 

-2aAN(a2jtN-m/p2)A(n)  (\0fa)  jplO^Np)^,  (aANrJ+l^,  (aANr)] 

+rK|n(aiNp)[l^_1  (aiNrJ+'^i  (aiNr)]} 

where  K^(ajtNp)  -  -  |  [l<m_ ,  (aiNpJ+K^ ,  (aXNp) 
and 

<»""■>  ‘  |  [V2(aJNr)+2'm<a*Nr>+W(a*Nr>] 


C-5 


(C-4b) 


(C-4c) 
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1 1 1  q*AN  <*  0 

rg2(u)-gp(-u)  ,  p  p 

|]m  - •-— = -  «  1  im  -J-(u-aAN)  (au-a  jtN+m/pd)  K  ((u+a!  A  N)p) 

u-0  L  u  J  u-«0  ul  m  ' 

•  eiu<x-*/a>A(">((W.)eb) 

*  (u+aUI N)  ((i«+a  | je|N)r)] 

+(-u-aiN)(™au-a2XM4fli/p2)Km((-u+aU|N)p)e",u^x"°r/a^A^((\+u/a)eb) 
•  Dm-i  ((-«*■•  M  |N)  r)  +  l|tH.1  ((-u+a  |i  |N)  r)  J 
-  2KJ3  |/N  |p)[lm_ ,  (a  |An|  rJ+l^j  (a|  IN|r)] 

•  -j2a2AN-m/p2-l(x-cr/a)aJeN(a2/N-m/p2)]A(n)  (\0fa) 

- !  9b^N (a2AN-m/p2) A  j (\0b)} 

-2a AN (a2AN-m/ p2) A (n)  (\eb){pl^<a|  AN|p)[lm-I(a|AN^  rj+l^,  (a|  AN|r)] 

+rKm(a|  AN  |p)[r  _ ,  (a  |  AN  | r)+  T _ ,  (a|  AN  |  r) ]]>  (C-4d) 


r>  /" 

C-c 
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APPENDIX  D 


Evaluation  of  the  Integrand  of  K  at  u  a  0 

9 


The  Integrand  of  K  ^m,n^ls 


9 


m 


(D-1! 


where  f(u)  -  (au-a2/N+  •—)  lm(lu-ajtNI  p)Km(  1  u-a  AN  I  r)e ' U  aVn*((X-  H)0b) 
q  *  jtN,  i  *  0,  ±  I ,  ±  2,  ... 

f(u)  *  (au  +  -Iy)lm(up)Km(ur)eiu(x"  a*A(n,((x-  £)0b) 

P 

f(-u)  -  (-auf  ^y)lm(up)Km(ur)e”lu(x“  a^A^n)((x+  -)9b) 

P 

Then 


llm  m[f (u)-f (-u)l  _  Ifm 
u-*0  U  Lh-0 


1  im 
u-«0 


m 


O  II 

2  m 


5  ln,<uP>Km(ur>{au[e'lJ'X"  j)6b) 

4  l>  A®  «X+X>] 

+  -^[elu<x*  |)A«f!)((x- i)eb) 

I 

-|u(x-  f>A<"'><(»  H)eb)]} 

a«im(up)Km(ur)[elu(x-  $  hM  i'.\-  j)3b) 

*  e-,u(*-  a>  A*"'*1**  X>] 
lin(up)Koi(ur)iVu(x'  f>  A(n)((\-  a)8b).a-|u(x‘  SW"> 


P 

e 


D-l 


X>] 
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a)  when  m  «  \  f  0 


HZ  KSaHkall  -  .<£)  A"K> 

2.  m  [e'“<x-  X>-«‘",(,‘  »)A(n)(u+  >b>] 

*  is  $  HS  1 - u - - 

-  a(A)"  A  (">  (X0b) 

+  TT  i  (r'  [' <x-  <A9b>  +  1  (x_  |)  A<n>  (Wb> 

.  .2  ,  ,*•  >b>  »A(!|(teXi 

7s  r  L  su  du  4-o 

-  $"  {[«K*-  f)4]  AW(Mb)  +  I  TT-^^KO  <»-'•> 


b)  when  m  -  \  ■  0,  the  Integrand  at  u  ■  0  Is 


llm  m[f(u)-f(-u)]  „ 
u-0 - z -  0 

In  fact,  the  entire  Integral  Is  zero. 


(D-lb) 


2)  q  -  JN  >  0 


ihO  ~UV*  U)'1  "  ilo  ;j{(au-a2£Wf  “j)  lm((aj2N-u)p)Km((aiN-u)r) 
.  e'u<x-  !>A(">((X-  f)9b) 

-  (-au-a22N+  -^)  lm((aiW-u)p)Kra(a«tnj)r) 

.  a-'“<x-f)  A(S)«0»b>} 

By  L’HospItal's  rule,  the  Integrand  at  u  -  0  is 


0-2 


li«  «!u(x“  5V  {[[»fl(x.-  |)(au-a2lHf  ^]lw((ajjNMi)p)Kw((alH-u)r) 
u-*0  P 

+  (au-a2|N+  T-  -f -p t pjK^CiajtM-u)  r) -r  I <a«N-u>  p>K^( (ajeH-u) r 


•  A(5)<a-X> 


+  (au-a2flt*  -!y)l||1((a«l-u)p)Km((atN-u)r) 


aA(")((x-T)9j 


9“ 


} 


r 

_e-lu(x- 2)m  JT  a_j  (x_  2.)  (-aU-a2jtN+  ay)  lm((ajgN+u)p)Km((ajjH+u)r) 

U.  p 

+  (-au-a2jtMf  ^)[pl^((ajlHfu)p)Km((ajlNfu)r)+rlm((ai»ku)p)K(||((aJtH4-a)r)] 


.  f)eb) 


+  (-au-a2jgN+  •2y)lm((alHfu)p)Km((aJ0*fu)r) 


3A(S)((x+  -)eb) 

du 


■} 


«=  2m  {[a+I(x-  |)(-aZXN+  y)]  lm(a  I  jfcNI  p)Km(a  J  IN  I  r)A(n)  (X0b) 

-  (-a2jtMf  -y)  A^n)  (X©b)[pl(J|^,»l -^WIpjKj^Cal  £NI  r)+rlm(al  iNI  p/ICj^al  INI  r)J 
P 

+  (-a2j&N+  -y)  lm(aljtNlp)Km(a6jlNlr)j^l  y  a|  ^  (XBj,]}  (d_2) 

P 

where 

i;(.l/Nlp)  -±[l„_,<aUHIp)  +  l^talANIp)] 

K^aljWlr)  *  -  |  [Km_,(ai iNlr)  +  K^falANIr)]  0>-2a) 

When  m  -  0  the  Integrand  Is  zero. 

3)  q  *  XN  <  0 

I  Im  m[f(u)  -  f (-u)  3  m 
u-tO  u 
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Hm  m  -f(au-aZxN  +  im((u-ajeN)p)Kn?((u-aJtN)r)elu^x“  a^A^((X-  -±)0b) 
u-*0  up' 

-  (-au-aZjfcN+  -^r)  lni((-u-aJtN)p)Km((-u-aJtN)r)e",U^X"  a^A^((\+  ~)eb)} 
P 

It  can  be  easily  shown  that 


I Im  m[f (u)  -  f (-u)] 
u-0  u 

2m  {[>!(x-  |)  (-aZjlN+  -^)]  Ijal  XNIpjKjallNI  r)  A(n)  (X0fa) 

p 

+  (-aZjtN+  *i|)A(n^(x9b)[pP(aliNlp)Kni(aljtNlr)+rlm(a|jtNlp)KiJi(aliNlr)] 

+  (-a2im-  -^)!fn(alXNIp)Km(alXNIr)[l  ^  (\9b)]  (0-3) 

P 

When  m  ■  0  the  Integrand  Is  zero. 


